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The status of possible new-physics signals in _B S!t j-meson mixing and decay is reviewed. In par- 
ticular, it is emphasized that the recent LHCb results, that find no evidence for a non-standard 
phase in B s -B 3 mixing, make a consistent explanation of the D0 data on the like-sign dimuon 
charge asymmetry notoriously difficult. In order to clarify the inconclusive experimental sit- 
uation, independent measurements of the semileptonic asymmetries are needed. 



1 Introduction 

The phenomenon of neutral l? s -meson mixing is encoded in the off-diagonal elements Mf 2 and 
Tf 2 of the mass and decay rate matrix. These two complex parameters can be fully determined 
by measuring the mass difference AM S = Mfj — M£, the CP-violating phase <f> s jj^ = — 2/3 s , the 
decay width difference Ar s = T S L — T S H , and the CP asymmetry a s j s in flavor-specific decays. 
The combined Tevatron and LHC determination of the mass difference readEE 

AM S = (17.73 ± 0.05) ps -1 , (1) 

and agrees well with the corresponding standard model (SM) predictiorP 

(AM S ) SM = (17.3 ±2.6) ps" 1 . (2) 

Here the quoted errors correspond to 68% confidence level (CL) ranges. 

The phase difference (pj/^p^ between the B s mixing and the b — > sec decay amplitude and the 
width difference Ar s can be simultaneously determined from an analysis of the flavor-tagged 
time-dependent decay B s — > J/ip<j). Such measurements have been performed by CDF, D0, 
and recently also by LHCb. Combining 1 fb -1 of B s — > J/yj(f> and B s — > J/ipfo data the LHCb 
collaboration findP 

#/W = (-°- 11±5 -°)° > Ar s = (0.116 ±0.019) ps" 1 . (3) 
The corresponding numbers in the SM arP 



J/i><t> 



ism = arg 



(V*V cb ) 2 



-2.1 ±0.1)°, (Ar s ) SM = (0.087 ±0.021) ps" 1 , (4) 



and again agree well with the observations with the results for Ar s differing by I.Oct. 

The CP asymmetry in flavor-specific decays a s j s can be extracted from a measurement of the 
like-sign dimuon charge asymmetry ^4 g L , which involves a sample that is almost evenly composed 
of Bd and B s mesons. Employing 9.0 fb -1 of data the D0 collaboration obtainP 

Ag L = (—7.87 ± 1.96) • lO -3 . (5) 

Utilizing the SM predictions for the individual flavor-specific CP asymmetrieP 

(o/Jsm = -(4.1 ±0.6) -10" 4 , (a^) S M = (1.9±0.3)-10- 5 , (6) 

one obtains 

04l)sm = (0.595 ± 0.022) (a d fs ) SM + (0.406 ± 0.022) (a s fs ) SM = (-2.4 ± 0.4) • 10" 4 . (7) 



Using the measured value of the CP asymmetry in flavor-specific decayP 



aJ fl = (-4.7±4.6)-l(r 3 , (8) 

and (|5j) in Q , one can also directly derive a value of a j s . One arrives at 

a s fs = (-1.3 ±0.8) • 1(T 2 , (9) 

The results §5§ and ([7]) correspond to a tension with a statistical significance of 3.9<r, while Q 
differs from (aj s )sM as given in (Tsj) by 1.5a only 



2 Model-Independent Analysis 

In view of the observed departures from the SM predictions, it is natural to ask what kind of 
new physics is able to simultaneously explain the measured values of AM S , <f>j/^, Ar s , and a^ s . 
This question can be addressed in a model-independent way by parametrizing the off-diagonal 
elements of the mass and decay rate matrix as follows 

M{ 2 = (Mf 2 ) SM + (Mf 2 ) NP = (Mf 2 )sM Rm e i4,M , 

(10) 

r i2 = ( r i2)sM + (rf 2 )NP = (r? 2 )sM Rr e^ r • 

In he presence of generic new physics, the l? s -meson observables of interest are then given to 
leading power in |rf 2 |/|Mf 2 | by 

AM S = (AM S ) SM Rm , <j>j/^ = 05/^)sm + Hi , 

cos (0| M + M - 0r) _ , Ar , D , x nu 

AF S = (Ar s )sM-Rr js ~ (Ar s ) SM Rr cos (<p M - <Pr) , (H) 

cos 4> S SM 

s , s \ sin (4> S SM + (p M -<fc) i s \ R r sin (4> M - 4>r) 

a fs = (O/JSM „ ^—5 ~ l«/JSM -5 Ti , 

3 3 Rm sm< ?sM Rm </>sm 

where the final results for Ar s and a s ^ s have been obtained by performing an expansion in 
0g M = arg (— (Mf 2 )sM/(rj 2 )sM)- Since 0g M = (0.22 ±0.06)°^ this is an excellent approximation. 

The four real parameters Rm,f and (f>M,r entering (10) can be constrained by confronting 
the observed values of AM S , <j>*j Ar s , and a s ^ s with their SM predictions. We begin our 
analysis by asking how well the SM hypothesis describes the data. Performing a global fit, 
we obtain x 2 = 3.5 corresponding to 0.7er (1.4<r) for 4 (2) degrees of freedom (dofs). For the 
best-fit solution, we find instead (Rm,4>m, Rr,<fir) = (1-02, 2°, 2.9, 65°). After marginalization 
the corresponding symmetrized 68% CL parameter ranges are 

R M = 1-04 ±0.16, <£ M = (-0.4±5.0)°, 
i?r = 3.4±1.8, r = (56 ± 22)° . 

Focusing on Rr and (pa, we see that a very good fit requires excessive corrections to Tf 2 . 

In order to further elucidate the latter point, we analyse two orthogonal hypothesis of new 
physics in B s -B s oscillations, namely a scenario with (M* 2 )np 7^ and (rf 2 )NP = and a 
scenario with (M* 2 )np = and (rf 2 )NP 7^ 0. The left (right) panel in Figure [T] shows the results 
of our global fit for scenario SI (S2) in the Rm cos 4>m~Rm sm(f>M (-Rr cos 4>y~Rt sinc^r) plane. 
From the left panel, we glean that in the scenario SI the regions of all individual constraints apart 
from a S j s overlap. Minimizing the x 2 function gives (R M , 4>m) = (1-03, 2.0°) and x 2 /dofs = 3.4/2 
corresponding to 1.3cr, which represents only a marginal improvement with respect to the SM 
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Figure 1: Left (Right): Constraints on Rm and <j>M (Rr and <f>r) m scenario SI (S2). For the individual 
constraints the colored areas are 68% CL regions (dofs = 1), while for the combined fit the red (light red) area is 
the 68% (95%) probability region (dofs = 2). The SM values (best-fit solutions) are marked by a dot (cross). 



hypothesisEE In consequence, the case for a non-zero non-standard contribution to M( 2 is rather 
weak. By inspection of the right panel, we see that in contrast to SI, in the scenario S2 a 
description of the data with a probability of better than 68% is possible. The best-fit point 
is located at (Rr,^>r) = (2.9,62°). In fact, the latter parameters lead to an almost perfect fit 
with x 2 /dofs = 0.2/2 corresponding to O.lcr. The data hence statistically favors the new-physics 
scenario S2 over the hypothesis SI. 



3 New Physics in b — > st+t 



The above findings suggest that one hypothetical explanation of the experimentally observed 
large negative values of a S f s (or equivalent ^4g L ) consists in postulating new physics in Tf 2 that 
changes the SM value by a factor of 3 or more. In the following we will study whether or not and 
to which extend such an speculative option is viable. While in principle any composite operator 
(sb)f, with / leading to an arbitrary flavor neutral final state of at least two fields and total 
mass below the _3<j-meson mass, can contribute to rf 2 , the field content of / is in practice very 
restricted, since B s — > f and — > X s f decays to most final states involving light particles are 
severely constrained. One notable exception is the subclass of B s - and _5rf-meson decays to a 
pair of tau leptons, as has been first pointed out a few years agoP 

The possibility of large b — > st + t~ contributions to Tf 2 , can be analyzed in a model- 
independent fashionr by adding 

C cS = 4 ^V*V tb Y,C l ^)Q l , (13) 

to the effective AB = 1 SM Lagrangian. Here [i denotes the renormalization scale and the Fermi 
constant Gp as well as the leading Cabibbo-Kobayashi-Maskawa (CKM) factor V t * s Vtb have been 
extracted as global prefactors. The index i runs over the complete set of dimension-six operators 
with flavor content (s6)(fr), namely (A, B = L,R) 

Qs,AB = (sP A b) (tP b t) , 

Qv,AB = (s-fPA b) (f lfl P B t) , (14) 
Q T , A = (sa^P A b)(fa flu P A T) , 




Figure 2: One-loop diagrams with a penguin insertion of a (s6)(rr) operator (orange squares) that contribute 
to the processes b — > s~/ (left), b s£ + £~ (middle), and b — > S77 (right). The tau loop in all graphs is closed. 



where Pl,r = (1 -F 7s) /2 project onto left- and right-handed chiral fields and a^ v = i [7^, 7"] /2. 



The ten operators entering (14) govern the purely leptonic B s — > t + t~ decay, the inclusive 
semi-leptonic B —> X s t + t~ decay, and its exclusive counterpart B + — > K + t + t~ , making these 
channels potentially powerful constraints. In practice, however, flavor-changing neutral current 
B s < i decays into final states involving taus are experimentally still largely unexplored territory, 
so that these direct constraints turn out to be not very strong. Explicitly, one obtains 

B(B s ^t + t-) < 3%, 
B{B^X s t+t-) < 2.5%, (15) 
B(B + -4 K+T+T-) < 3.3 • 1(T 3 . 

Here the first limit deriveP from comparing the SM prediction t b Jt b<1 - 1 £ [-0.4, 0.0]%E1 with 
the corresponding experimental result t Bs J T Bd — 1 = (0.4±1.9)%, while the second (crude) bound 
follows from estimating the possible contamination of the exclusive and inclusive semileptonic 
decay samples by B — > X s t + t~ events. Limits on B(B S — > t + t~) and B(B — > X s t + t~) of 
strength similar to those given in ( |l5[ ) also follow from charm counting and/or LEP searches 
for B decays with large missing energy^ The final number corresponds to the 90% CL upper 
limit on the branching ratio of B + — > K + t + t~ as measured by BaBarP^l 

Further constraints on the Wilson coefficients of the (sb)(fr) operators arise indirectly from 
the experimentally available information on the b — > S7, b — > s£ + £~ {I = e,fi), and b — > 577. 



transitions, because some of the effective operators introduced in (14) mix into the electromag- 
netic dipole operators (J7,A and the vector-like semileptonic operators Qq,a- The relevant Feyn- 
man diagrams are shown in Figure [2J An explicit calculatioiP shows that the operators Qs,AB 
mix neither into Qt^a nor Qg t A, while Qv,AB (Qt,a) mixes only into Qg,^ {Qt,a)- As a result 
of the particular mixing pattern, the stringent constraints from the radiative decay B — > X s ^f 
rule out large contributions to T\ 2 only if they arise from the tensor operators Qt.a- Similarly, 
the rare decays B — > X s £ + £~ and B — > K^t^t' primarily limit contributions stemming from 
the vector operators Qv,AB- In contrast to B — > X s ^y, all (s6)(fr) operators contribute to the 
double-radiative B s — )■ 77 decay at the one- loop level. A detailed study^ shows however that the 
limits following from b -4 S77 are in practice not competitive with the bounds obtained from 
the other tree- and loop-level mediated i3 Sj rf-meson decays. 

The model-independent 90% CL limits on the magnitudes of the Wilson coefficients are 
summarized in Table [TJ We see that in the case of the (s6)(rr) scalar and vector operators 
the allowed effects can reach almost 0(1). These Wilson coefficients can hence be similar 
in size to that of the color-singlet current-current operator Q2, which provides the dominant 
contribution to Tf 2 in the SM. In consequence, the corresponding suppression scale A is quite low, 
around (1 — 2) TeV. Possible new-physics contributions to the (s6)(rr) tensor operators are more 
severely constrained, because they lead to a contamination of B — > X s ^f. It is also interesting 
to ask which impact future improved extractions of B(B S — > r + r _ ), B{B — > X s t + t~), and 
B(B + — > K + t + t~) will have on the limits on the Wilson coefficients C^m^). Such a comparison 



Operator 


Bound on d(mb) 


Bound on A 


Observable 


(sP A b)(fP B T) 


0.5 


2.0 TeV 


B s r+r" 


(s^P A b){f^P B T) 


0.8 


1.0 TeV 


B+ -> K+r+r- 


(sa^P L b)(fa^P L T) 


0.06 


3.2 TeV 


B -> X s7 


{sa^P R b)(fa^P R T) 


0.09 


2.8 TeV 





Table 1: Model-independent limits on the Wilson coefficients of (s6)(rr) operators. The second (third) 
column shows the limit on Ci(mb) (the bound on the suppression scale A assuming a coupling strength of 1). 



is provided in Figure[3j From the left panel one concludes that B(B + — > K + t + t~) < 1.3 • 10~ 3 , 
corresponding to an improvement of the present upper limit by a factor of 2.5, would allow to set 
a bound on Cs,A(mb) that is as good as the one that follows at present already from B s — > t + t~ . 
As can seen from the right panel, in the case of Cy; AB (?™&) ; limits of B(B S — >■ t + t~) < 1.7% and 
B(B —> X s t + t~) < 0.7% are needed to be competitive with B + — > K + t + t~ . The quoted limits 
correspond to improvements of our current knowledge of the relevant S-meson branching ratios 
by a factor of 1.8 and 3.6, respectively. Finally, for what concerns CT,A( m b)> even improvements 
of the direct constraints by a factor of more than 10 are not sufficient to beat the indirect 
constraint arising from B — > X s ^. 



Effects of (sb)(rr) Operators in T s 



12 



The off-diagonal element of the decay-width matrix is related via the optical theorem to the ab- 
sorptive part of the forward-scattering amplitude which converts a B s into a B s meson. Working 
to leading order in the strong coupling constant and Aqcd/^6, the contributions from the com- 
plete set of operators (14) to T\ 2 is found by computing the matrix elements of the (Qi,Qj) 



double insertions between quark states. Such a calculation^ leads to the results 

(Rr)s,AB < 1 + (0.4 ± 0.1) \C s ,AB(m b )\ 2 , 

(Rr)v,AB < 1 + (0.4 ± 0.1) \C v ,AB(m b )\ 2 , (16) 

(Rt)t,a < 1 + (0.9 ± 0.2) \C T ,A( m b)\ 2 , 

where the quoted uncertainties are due to the error on (Ar s )sM as given in Q. Employing now 
the 90% CL bounds on the low-energy Wilson coefficients given in Table [TJ it follows that 

(Rr)s,AB < 1-15 , (-Rr)v,AB<1.35, (Rv)t,l < 1-004, (R v )t,r < 1-008 . (17) 

These numbers imply that (s6)(rr) operators of scalar (vector) type can lead to enhancements 
of |rf 2 | over its SM value by 15% (35%) without violating any existing constraint. In contrast, 
contributions from tensor operators can alter |rf 2 | by at most 0.8%. These numbers should be 



compared to the best-fit solution for Rr as given in (12). From the comparison it immediately 
becomes apparent that absorptive new physics in Tf 2 in form of (sb)(fr) operators cannot provide 
a satisfactory explanation of the anomaly in the dimuon charge asymmetry data observed by 
the D0 collaboration. This is a model-independent conclusion that can be shown to hold in 
explicit models of new physics with modification of the b — > st + t~ channel such as leptoquark 
scenarios or Z' models'? 



5 Conclusions and Outlook 

Motivated by the observation that the anomalously large dimuon charge asymmetry measured 
by the D0 collaboration, can be fully explained only if new physics contributes to the absorptive 
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Figure 3: Left (Right): Bound on Cs,A{iTit) (Cv,AB(m,b)) from improved determinations of B(B S — > t + t~) (red), 
B(B — ¥ X s t + t~) (green), and B(B + — > K + t + t~) (blue). The dashed lines indicate the presently best bound. 



part of the B s ^~B s ,d mixing amplitudes, we have presented a model-independent study of the 
contributions to Tf 2 arising from the complete set of dimension-six (s6)(fr) operators. Taking 
into account the direct bounds from B s — > t + t~ , B — > X s t + t~ , and B + — > K + t + t~ as well 
as the indirect constraints from b — > 57, b — > s£ + £~, and b — > S77, we have demonstrated 
that only the Wilson coefficients of the tensor operators are severely constrained by data, while 
those of the scalar and vector operators can be sizable and almost reach the size of the Wilson 
coefficient of the leading current-current SM operator. It follows that the presence of a single 
(s6)(rr) operator can lead to an enhancement of Tf 2 of at most 35% compared to its SM value. 
Since a resolution of the tension in the i3 s -meson sector would require the effects to be of the 
order of 300% (or larger), the allowed shifts are by far too small to provide an satisfactory 
explanation of the issue. We emphasize that after minor modifications, our general results can 
be applied to other dimension-six operators involving quarks and leptons. For example, as a 
result of the 90% CL limit B(B + -)■ 2T+Y ± ^ T ) < 7.7 • 10~ 5 ^ the direct bounds on the Wilson 
coefficients of the set of (sb)(ffi) operators turn out to be roughly a factor of 7 stronger than 
those in the (s6)(rr) case. Possible effects of (sb)(f[i) operators are therefore generically too 
small to lead to a notable improvement of the tension present in the current 5 s ^-meson data. 
Similarly, a contribution from (d6)(rr) operators to Tf 2 large enough to explain the ^4g L data 
is excluded by the 90% CL bound B(B -)■ t + t~) < 4.1 • 10~ 3 P Naively, also (bd)(cc) operators 
are heavily constrained (meaning that their Wilson coefficients should be smaller than those of 
the QCD/electroweak penguins in the SM) by the plethora of exclusive B decays. A dedicated 
analysis of the latter class of contributions is however not available in the literature. 

Our model-independent study of non-standard effects in Tf 2 can readily be applied to explicit 
SM extensions involving leptoquarks or Z' bosons. In fact, the pattern of deviations found in 
these scenarios resembles the one of all new-physics model with real tnp = (Mf 2 )Np/(rf 2 )NP, for 
which it can be shown that the measurement of AA1 S generically puts stringent constraints on 
both Ar s and a s j s . These bounds turn out to be weakest if the ratio tnp is positive and as small 
as possible. Since on dimensional grounds tnp scales as the square of the new-physics scale, this 
general observation implies that SM extensions that aim at a good description of the Tevatron 
data should have new degrees of freedom below the electroweak scale and/or be equipped with a 
mechanism that renders the contribution to Mf 2 small. Furthermore, models in which (M{ 2 )np 
is generated beyond Born level seem more promising, since in such a case r^p is suppressed by 
a loop factor with respect to the case where (Mf 2 )NP arises already at tree level. 



The above discussion implies that a full explanation of the observed discrepancies is not even 
possible for the most general case (Mf 2 )NP 7^ and (rf 2 )NP 7^ 0- Numerically, one finds that 
the addition of a single (s6)(rr) vector operator giving (Rt)v,ab = 1-35 on top of dispersive 
new physics with (M^ 2 )np 7^ 0, can only improve the quality of the fit to the latest set of 
measurements to \ 2 = 1-4 compared to x 2 = 3.5 within the SM. This might indicate that 
the high central value of observed by the D0 collaboration is (partly) due to a statistical 
fluctuation. Future improvements in the measurement of the CP phase (pj/^^ and, in particular, 

a first determination of the difference a s j s — a^ s between the B s and semileptonic asymmetries 
by LHCb, are soon expected to shed light on this issue, and are of utmost importance in order 
to answer whether or not there is new physics hiding in the -B^-meson sector. 
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